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Introduction 


With the evolution of advanced composites, the feasibility of designing bearingless rotor systems 
for high speed, demanding maneuver envelopes, and high aircraft gross weights has become a 
reality. These systems eliminate the need for hinges and heavily loaded bearings by incorporating 
a composite flexbeam structure which accommodates flapping, lead-lag, and feathering motions by 
bending and twisting while reacting full blade centrifugal force, Figure (1). The flight characteristics 
of a bearingless rotor system are largely dependent on hub design, and the principal element in this 
type of system is the composite flexbeam. As in any hub design, trade off studies must be performed 
in order to optimize performance, dynamics (stability), handling qualities, and stresses. However, 
since the flexbeam structure is the primary component which will determine the balance of these 
characteristics, its design and fabrication are not straightforward. Some of the considerations 
which must be addressed are as follows: 

1. Flap-Lag-Torsion deformations will be accommodated through the flexbeam. 

2. Effective flapping hinge offset has to be properly controlled for a balance between maneuver- 
ability and dynamic vibration. 

3. Hub size must be kept at a minimum in order to reduce weight and hub drag. 

4. Optimum tailoring of the pitchcase, snubber/damper and inplane flexbeam deformation must 
be obtained in order to maximize inplane damping. 

5. The flexbeam must be able to endure peak loading from high g maneuvers as well as endurance 
flight loads. 

6. Flexbeam design criteria are influenced by rotor shaft /mast/hub impedance characteristics. 
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Figure (1) 



Model Development 


At McDonnell Douglas Helicopter Company, the previous considerations have been integrated into 
a mathmatical procedure for design and optimization of advanced composite bearingless rotor 
systems. Because of the highly coupled structural and damping requirements of the flexbeam, 
it was critical to include a representation of the pitchcase and snubber/damper, as well as the 
flexbeam, in the math model. Various alternatives can be proposed for formulation of such a 
problem: 

1. choice of flexbeam cross section 

2. choice of composite material type and configuration 

3. choice of hub attachment configuration 

4. choice of objective function and constraints for optimization routine 

(a) minimization of peak stresses due to flap, lag and torsion with a requirement on damper 
motion and hence rotor inplane damping 

(b) maximization of rotor damping with upper bound contraints on stresses 


Flexbeam Cross Section 

In an effort to simplify the beam design and fabrication, a comprehensive review of candidate 
cross sections was conducted under a company IRAD program . Geometries studied ranged from 
rectangular, cruciform and multiple “H” closed sections to multiple element open sections and 
elastomericly connected rectangular cross sections. 

Based on this review, a simple rectangular section was selected as the best overall configuration 
for the advanced flexbeam. In addition to the efficient load carrying capability, the relatively 
straightforward rectangular configuration offered advantages of simplified fabrication resulting in 
lower production costs, superior quality control, and simplified inspectability. 


Composite Material Selection 

The composite material selected for the advanced flexbeam was S-2 fiberglass impregnated with a 
350 deg F cured “toughened resin” epoxy system. The resin system was selected based on the results 
of an extensive study of candidate systems conducted under other McDonnell Douglas Helicopter 
Company IRAD programs. 

The S-2 fiberglass provides improved compressive stress margins and reduces long term enviromen- 
tal degradation concerns associated with the previously selected Kevlar material. 


* Internal Research and Development 
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Hub Attachment Configuration 


Flexbeam design criteria are influenced by rotor shaft/ mast/hub impedance characteristics. A 
relatively soft rotor mast mount can significantly modify a flexbeam’s effective flapping hinge offset 
and bending requirements. An inhouse analysis has indicated that a softer mount can result in an 
increased level of airframe vibration and the need for vibration suppression systems. 

The Helicopter Advanced Rotor Program (HARP) flight test program demonstrated low vibra- 
tion characteristics for the initial composite flexbeams using the McDonnell Douglas Helicopter 
Company rigid mast system with no active or passive vibration control equipment required. The 
advanced beam design approach was also based on the proven rigid mast configuration. 


Objective Function and Constraints 

Significant progress has been made both by industry and the goverment in the study of flexbeam 
design over the past ten years. McDonnell Douglas Helicopter Company has kept abreast of this 
work and has incorporated the lessons learned both by ourselves and other investigators into the 
present math model. 

The primary goal in flexbeam sizing is to achieve the required rotor blade motions without exceeding 
the allowable stresses in the composite material. Simultaneous goals are to minimize flexbeam 
length, to achieve blade flap and lead-lag hinge offsets which provide suitable dynamics properties, 
and to obtain sufficient lead-lag motion at the snubber/damper for satisfactory dynamic response 
to cyclic loads. The criteria for blade motions include the maximum blade cyclic flap, lead-lag, and 
feathering motions seen in the flight envelope (the one hour condition), as well as endurance limit, 
static droop and start-up conditions. 

Rotor blade flapping, lead-lag and feathering motions acting in combination with the blade and 
pitchcase centrifugal force cause flexbeam deformations. The flexbeam must be so designed as to 
allow these blade motions without exceeding the allowable stresses of the fiber composite of which 
it is constructed. Because blade motions have both steady and cyclic components, fatigue loads 
prevail and the fatigue strength of the composite material is critical for design. Therefore, the one 
hour and endurance limit flap, lead-lag, and torsional motions are applied to the flexbeam. The 
cyclic stresses which the one hour motions cause should not exceed the corresponding one hour 
fatigue allowables of the composite material. Similarly, cyclic stresses due to endurance limit loads 
should not exceed endurance limit fatigue allowables. In addition, criteria for the blade motions 
include static droop and start-up torque conditions. These goals, taken together, are used in sizing 
the flexbeam in lieu of a complete fatigue life analysis covering the full operating spectrum of the 
helicopter. 

The objective function selected for the advanced beam was minimization of stresses due to a critical 
combination of flap, lead-lag and torsional deformation. In addition, a displacement constraint is 
imposed for adequate damper motion to satisfy the dynamic requirements. 
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Theoretical Background 


The conventional finite element method is not applicable in the case of an elastic structure under- 
going high angular motions. This is of interest because numerous structural configurations such 
as spinning helicopter blades, satellites, rotating flexbeams, shafts and linkages fall into this cate- 
gory. The analysis of these rotating structures differs from that of stationary structures due to the 
complexity of the accelerations which act throughout the system. In addition to the accelerations 
resulting from elastic structural deformations, contributions due to Coriolis and centripetal accel- 
eration may be of significance. Also, the stiffness characteristics of the structure may be modified 
by the internal loads induced by the centrifugal forces. The finite element formulation given in this 
section includes all these effects in a uniform general formulation. 

The variational principle useful for problems in dynamics is Hamilton’s principle. The functional 
for this principle is the Lagrangian, L, defined as 


L = T — U — V (1) 

where T is the kinetic energy, U is the strain energy, and V is the potential of the applied loads. 
The strain energy density may be expressed as 

dU = \e T adV = \e T CedV (2) 

Z Z 

where C is the material elastic matrix. The potential of the applied loads, V, may be expressed as 
the sum of the potential of the body forces and surface tractions as 

V = - j v (X -u)dV - j s {T -u)dS (3) 

where X is the body force vector, u is the displacement vector, and T is the vector of surface 
tractions. The second integral is over the surface of the body on which the surface tractions are 
prescribed. 

The total kinetic energy density is defined as 


dT = - P V T VdV (4) 

z 

where p is the density of the material and V is the absolute velocity of a particle in the element as 
given by 


V — Vrb + u + Cju 


( 5 ) 
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Here, V RB is the velocity of the particle due to frame motion, u is the vector of elastic deformations, 
and u> is the angular velocity of the element coordinate system. Substitution of equation (5) into 
(4) gives 


dT = \pV* B V RB dV + l 2 pu T udV 
+ Jf(u r wu — u T u)u — u T CjCju)dV 

Li 

+ + u t V rb + V RB wu - u T u)V RB )dV 


(6) 


The first term in equation (6) is the kinetic energy due purely to frame translation. This term will 
be deleted since frame motion is assumed to be entirely prescribed. The second term represents 
the kinetic energy due to elastic deformations only. The third term is the energy due to frame 
rotation, and the fourth term gives the forces due to coupling between frame translation, elastic 
deformations and frame rotation. 

Using equations (2), (3) and (6), the functional in equation (1) may be expressed in the following 
form: 

L — i J p (u T u + u t Cju — u T d)u — u T iiwu) dV 

-- [ e T CedV 
2 Jv 

+ \!v p ^ Vrb ^ + ^ TVrb + v nB^ u ~ u T CjV RB )dV 

+ J u T XdV + J u T TdS (7) 


The statement of Hamilton’s principle is as follows: Among all possible time histories of displace- 
ment configurations which satisfy compatibility and the constraints or kinematic boundary condi- 
tions and which also satisfy conditions at times ti and t 2 , the history which is the actual solution 
makes the Lagrangian functional a minimum. This principle may be expressed mathematically as 


ft 2 

6 Ldt = 0 (8) 

Jh 


Application to the Finite Element Method 

Most problems are too complex to use energy principles to obtain the exact solution directly. 
The usual technique is to guess, a priori, a trial family of solutions for the unknown quantity in 
order to construct a variational functional. We therefore choose an approximate pattern for the 
trial solution and apply the energy principles to the approximate functional to obtain the equations 
necessary to find the approximate solution. In the theory of the finite element method for structural 
analysis, once the proper variational principle has been selected for the given problem, we express 
the functional involved in terms of approximate assumed displacement functions which satisfy the 
geometric boundary conditions. We then minimize the approximate functional using equation (8) 
to obtain a set of governing equations. 
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General Displacement Models 


The basic philosophy of the finite element method is piecewise approximation. That is, we approx- 
imate a solution to a complicated problem by subdividing the region of interest and representing 
the solution within each subdivision by a relatively simple function. In the displacement method of 
structural analysis, the structure or body is divided into finite elements. Then, simple functions are 
chosen to approximate the displacements within each element. These functions are called displace- 
ment models, displacement functions, displacement fields, or displacement patterns. A polynomial 
is the most common form of displacement model used for two principal reasons. First, it is easy to 
handle the mathematics of polynomials in formulating the desired equations for various elements 
and in performing digital computation. In particular, the use of polynomials permits us to dif- 
ferentiate and integrate with relative ease. Secondly, a polynomial of arbitrary order permits a 
recognizable approximation to the true solution; however, for practical purposes we are limited to 
one of finite order. By truncating an infinite polynomial at different orders, we clearly vary the 
degree of approximation. 

An exact solution for the displacement u(x ) is then approximated by various degree polynomials 
of the general form 


u{x) = ai + a 2 x + a 3 x 2 + 1- a n+ ix n (9) 

The greater the number of terms included in the approximation, the more closely the exact so- 
lution is represented. In equation (9), the coefficients of the polynomial, the a’s, are known as 
generalized coordinates or generalized displacement amplitudes. The number of terms retained in 
the polynomial determines the shape of the displacement model, whereas the magnitudes of the 
generalized coordinates govern the amplitude. These amplitudes are called generalized because 
they are not necessarily identified with the physical displacements of the element on a one-to-one 
basis; rather, they are linear combinations of some of the nodal displacements and perhaps of some 
of the derivatives of displacements at the nodes as well. The generalized coordinates represent the 
minimum number of parameters necessary to specify the polynomial amplitude. 

Equation (9) can be expressed in vector form as 

u(x) — <f> T a (10) 


where 

(f) = { lxx 2 x 3 • • • x n } T , 


and 


a = {a 1 a 2 «3 ' 1 • « n +i} T ' 


The degrees of freedom can thus be related to the generalized coordinate system by employing the 
displacement model. We can evaluate the generalized displacements at the nodes by substituting 
the nodal coordinates into the model. For example, using a model of the form given by equation 
(10), we may write, for a single node point 


u(x) = <j>a 


( 11 ) 
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where u is the vector of degrees of freedom for a single node; or, for the entire element, 


( u(node 1) 


<f>(node 1) 

1 u(node 2) 

X — \ 

► = 

4>{node 2) 

( u(nodeN n ) t 


4>{nodeN n ) 


( 12 ) 


where N n is the total number of nodes for the element being considered, x is the vector of nodal 
degrees of freedom, and the notation in parentheses indicates that the dependent variables are 
assigned their values at the particular node. We may invert equation (12) to get 

a = A~ x x (13) 

where A~ l is a displacement transformation matrix. Note that A is a square matrix, hence, the 
total number of generalized coordinates equals the total number of joint and internal degrees of 
freedom. Equation (13) may then be used to eliminate the generalized coordinates to obtain 

u = <f>A~ l x = Nx (14) 


which expresses the displacements u at any point within the element in terms of the displacements 
of the nodes x. 

Since the derivations of the element matrices are performed in the a coordinate system, we fur- 
ther develop the necessary relationships in terms of the generalized a coordinates rather than the 
physical x coordinates. 

The strains are expressed in terms of some combination of the derivatives of the displacements u. 
Since the generalized coordinates a are not functions of the spatial coordinates, these derivatives 
must be performed in terms of the matrix 4>. If s is the vector of the relevant strain components 
at an arbitrary point within the finite element, we use the strain-displacement equations and the 
displacement model to write 


e = B a oc (15) 

It is important to note that for geometrically nonlinear problems e and B a are functions of the 
independent space coordinates as well as functions of the generalized coordinates. Decomposing 
these quantities into linear and nonlinear parts, we have 


e — ej + e„i 


and therefore, 


B a = B ai + B 


a„t 


(16) 


(17) 
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If a is the vector of stresses corresponding to the strains e, we may use an appropriate matrix form 
of the stress strain equations and equation (15) to write the element stresses as 


a = CB a a (18) 

where C is the matrix of material constants given by equation (3). 

Variational Formulation of Element Matrices 

As discussed above, the method used here for calculating the element matrices and load vectors is 
the application of Hamilton’s principle. Therefore, we develop the Lagrangian given in equation (7) 
in terms of the a generalized coordinates and apply equation (8) to obtain the element equations 
in terms of the generalized coordinates and shape functions. 

In order to formulate the Lagrangian in terms of the generalized coordinates, we insert equations 
(11), (15) and (18) into equation (7) to obtain 

L=~ f a T B T a CB a ctdV 
2 Jv 



— ^ J p {cc T <j> T 4>a + a T B 2 a — a T B 2 a — a. T Biaj dV 

~\fv P + + V£ B u<t>a - <x t 4> t CjV R b) dV 

- j a T d> T XdV - j a T <t> T TdS 

(19) 

where 


B\ = <f> T Cjh)<f> 

(20) 

and 


B 2 — <f> T w<f> 

(21) 


Applying the variational principle given by equation (8), we obtain 

J* 2 (Sc? f v {B T a CB a + B a )dVa - Sa T J v p<j> T <j>dV ex') dt 
+ J* 2 (~6a T j y P B 2 dVa + Sa T pB 2 dVc dt 

+ (do? J pB x dV cx) dt 

— Jt (6a T p<j> T V RB dV — Sa T p(j> T u>V RB dv') dt 
- J* 2 ( 6a T j v (j> T XdV + 6a T jf 4> T Tds) dt = 0 (22) 
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where the matrix B a arises from the variation of the nonlinear part of the matrix B a which is a 
function of the generalized coordinates a. The t th row of the matrix B a is given by 


B a; = <x T Bl.C 


,dB anl . 


(23) 


a, ~ dcti 

for t = 1,2, • • • , 12. Integration of the terms in equation(22) by parts with respect to time gives 

j*' 6a T p<f> T <f>d.V & + jT 2 pB 2 dVcrj dt 

+ J*' Sa T ( f y {BlCB a + B a + pB 2 + pBjdVa') dt 

+ 6a T {Jy^P^ T V rb + p<f> T &VjiB)dV ^ dt 

- J*' Sa T <f> T XdV - j s <f> T TdS^j dt = 0 (24) 

Since the variations of the generalized displacements, 6a, are arbitrary, the sum of the expressions 
in parentheses must vanish. Therefore we obtain the equations of motion for the element as 


M“d + C a oc + K a a = F a (25) 

where M a is the consistent mass matrix defined by 

M a = [ p<j> T <j>dV (26) 

Jv 

C a is the damping matrix due to frame motion as given by 

C a = 2 J v pB 2 dV (27) 

and K a is the element stiffness matrix in terms of a coordinates as given by 

K a = J v {BlCB a + B a )dV + J v p(Bi + B 2 )dV (28) 

The first integral in the above equation is the stiffness due to strain energy and will be referred to 
as if“. The second integral is the stiffness due to frame motion to be referred to as K^. Recall 
from equation (17) that B a is comprised of linear and nonlinear parts. Then the stiffness due to 
strain energy, K“, may be decomposed into linear and nonlinear parts to give 

K = l BlCB a ,dV (29) 
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(30) 


Xt, = l BlCB.„dV + f y Bl,CB a ,dV 

+ 1 B„JCB..,dV + j y B.dV 

The stiffness matrix resulting from the first three terms in equation (30) is known as the large 
displacement matrix. The stiffness matrix arising from B a is dependent on the stress level and is 
known as the initial stress matrix, Ref. (l). 

Finally, the load vector, is defined by 

F a = J v <f> T XdV + j <t> T TdS 

- f v {p<t> T V RB + p<t> T CjV RB )dV (31) 

Displacement Models for a Beam Element 

Having formulated the expressions for the element matrices and loads in terms of the general shape 
functions, the next step is to choose the displacement models that will approximate the solution 
for the problem at hand. A flexbeam is a cantilever beam under tension. For a beam element, 
axial, lateral and torsional deformations are of interest. To approximate the deformation field in 
a flexbeam, the polynomial shape function of equation (9) is truncated to obtain the approximate 
representation of the true axial, bending and torsional displacements. 

Linear and cubic shape functions are used to express the displacements of the beam in terms of the 
a generalized coordinates. For axial deformations due to axial forces, the following linear shape 
function is employed: 

«<» = <*! + a^x (32) 

This shape function satisfies rigid body motion, constant strain states conditions and the compat- 
ibility conditions. 

For deformations due to bending about the y and z axes, the following cubic shape functions are 
used: 

^ = a 3 + a 4 x + a B x 2 + a 6 x 3 (33) 

tv b = a 7 + a 8 x + a 9 x 2 + a i0 x 3 (34) 

where is the displacement due to bending about the z axis and is the bending displacement 
about the y axis. These shape functions are also compatible and complete. The compatibility 
conditions at the ends of the element are met on displacements as well as the slopes. 
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For deformations due to torsion, the following linear shape function is used 


0 = an + ai 2 x (35) 

With the above shape functions, and the general expressions for the element matrices given by equa- 
tions (26) to (28) and (31), the element matrices can be easily obtained via a symbolic manipulation 
program such as SMP. 

Having the element matrices constructed, the assembly and the solution of the equations of motion 
are a routine matter. 


Simplification of Equations of Motion 

Element matrices given by equations (26) to (30) can be further simplified for the flexbeam appli- 
cation. We start with equation (28), the second integral, which represents the softening effect of 
the frame rotation on the element stiffness matrix. This term is usually small compared to the first 
integral if the rotor angular velocity is not close to one of the structural natural frequencies. For 
flexbeam design, it is small if the rotor angular velocity is not close to the first lead-lag frequency of 
the nonrotating flexbeam. Neglecting this term will make the model stiffer than the real flexbeam, 
and hence will make the prediction of stresses more conservative. 

The second set of terms which can be neglected in the flexbeam design is the first three integrals in 
equation (30). These integrals represent the effect of large rotations on the equilibrium equations 
of the flexbeam. For angles smaller than 15° this term is much smaller than the structural stiffness 
,ff“ and the centrifugal stiffness B a . In a flexbeam design this term can also be neglected, resulting 
in more a conservative design. 

Finally, the calculation of the centrifugal stiffness matrix, B a , can be simplified by noting that most 
of the CF force is caused by the weight of the blade and the pitchcase. Since the weight of the 
flexbeam is negligible compared to the weight of the blade and the pitchcase, we can assume that 
the CF force remains constant throughout the length of the flexbeam. The magnitude of this force 
can then be calculated separately and used as input to the program. This will eliminate the extra 
calculations otherwise needed to compute this force and hence results in saving computational time. 
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Equivalent Flexbeam Composite Properties 


The composite properties of the flexbeam in the principal material directions can be obtained using 
the Halpin-Tsai equations of Ref.(2) as follows: 


Ei = $E f + (1 - $)E m 

(36) 

_ (1 + 2f) 

E 2 (i - <r) Em 

(37) 

£ _ G m ( Gf + G m + ( Gf — G m )f] 

Gf + G m — (Gf — G m )$ 

(38) 

Vl2 = + Vm( 1 - f) 

(39) 

where 


Ei — composite Young’s modulus in fiber direction 
E 2 = composite Young’s modulus pependicular to the fiber direction 
(7 12 = composite shear modulus 
= composite Poisson’s ratio 

Ef , E m = Young’s moduli of fiber and matrix, respectively 
G f , G m = shear moduli of fiber and matrix, respectively 
i/ f , u m — Poissons’s ratio of fiber and matrix, respectively 



£ = fiber volume fraction 


Rigid Constraint Element 

For bearingless rotor flexbeam design, a pitchcase model is required. One might attempt to use a 
stiff beam element for the pitchcase model; however, a stiff beam will result in an ill-conditioned 
stiffness matrix. Since a kinematic constraint involves a relationship between degrees of freedom, 
the correct approach is that of employing a multipoint constraint. 

For large angular motion constraints between the degrees of freedom of an element, ideally a 
nonlinear constraint relationship is required. However, in flexbeam design, the angular rotations of 
flexbeam and pitchcase are small, and hence a linear constraint is sufficient. 
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To derive the rigid constraint element which will be used to represent the pitchcase, consider two 
node points A and B , which are connected by a rigid link. The rotational degrees of freedom of 
nodes A and B are 0 A and 0 B , respectively. Similarly, the displacement degrees of freedom are u A 
and u B . Assuming a rigid member between the two nodes, the displacements and rotations of node 
B can then be expressed in terms of the displacements and rotations of node A as 

u B = u A + 9a x r BA (40) 


and 


6 a — Ob 


(41) 


where r BA is a vector from node A to node B the magnitude of which is the length of the rigid 
member. Writing the vector cross product in matrix form, we obtain 


•uf 


•i 

0 

0 

0 

L t 

L v 1 


'Uf 

u 2 


0 

1 

0 

-L. 

0 

L, 


U2 

u s 


0 

0 

1 

L v 

-L x 

0 


U 3 

0i 


0 

0 

0 

1 

0 

0 


*1 
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0 

0 

0 

1 

0 


o 2 

- 03 ■ 

B 

.0 

0 

0 

0 

0 

1 


- #3 - 


(42) 


This, then, is the constraint equation relating the degrees of freedom of the outboard end of the 
fl exbeam to the degrees of freedom of the snubber /damper element. 


HUBFLEX Mathematical Model 

Using the mathematical formulations presented earlier, a finite element computer program called 
IIUBFLEX was developed. This computer program is tailored specifically for the analysis of 
flexbeams with rectangular cross sections. 

IIUBFLEX is a finite element analysis model for a cantilevered beam under centrifugal force. The 
model permits rapid calculation of spanwise load and stress distributions for a specific geometry 
flexbeam with specified material properties, centrifugal force, and modal deflections. The analysis 
has been validated against two- and three-dimensional NASTRAN models of equivalent beams, and 
has shown excellent correlation. Spanwise stress distributions can be calculated with HUBFLEX 
in less than 10 seconds CPU time on a Digital Equipment Corporation VAX 11/785. Substantial 
economies in the study of new designs are thus realized. 

Analytically, HUBFLEX treats the flexbeam as fixed at the inboard end of the rotor mast. At 
the outboard end, radial tension, bending moments, shears, and torque are applied. The beam 
is a statically indeterminate structure in that the elastic/cross sectional properties influence the 
load distribution. When any section of the beam is modified, a complete iteration is required to 
determine the new load and stress distribution along the beam. 
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The IIUBFLEX model utilizes the stiffness matrix method presented earlier to solve for flexbeam 
load distributions. Tension beam effects are modeled using centrifugal stiffening terms in each 
element, B a in equation (28). The effect of the pitchcase redundant load path is modeled by 
a rigid element, equation (42), attached to the flexbeam at the outboard end and attached to 
the flexbeam at the inboard end by a snubber/damper of given stiffnesses in the flapwise and 
chordwise directions. In the chordwise direction, the spring stiffness corresponds to the stiffness of 
the damper pads. In flap, the spring rate is equal to the stiffness of the snubber corrected for the 
pitchcase flapwise flexibility. The pitchcase is rigidly attached to the flexbeam at its outboard end. 
Inputs to the program include section element definition, flexbeam geometry, material properties, 
applied loads, and the desired deflections/rotations at the blade attachment point. The input forces 
and moments represent those applied to the flexbeam by the blade spanwise moment and shear 
distributions. The shear to moment ratios at the blade attachment have a powerful influence on 
the flexbeam deflected shape and are determined separately in an iterative procedure. They are 
continually updated as the configuration evolves. The spanwise load and stress distribution results 
include the effect of combined loadings associated with the relative phasing of the flap, chord, and 
torsion motions. As with the shear to moment ratios, the phasing has a powerful influence on the 
final results especially for the critical shear loads due to combined normal bending and feathering. 
The phasing relationships were defined using an aeroelastic analysis (DART) and were verified by 
the HARP flight test results. 
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HUBFLEX/ADS Computer Program 


The HUBFLEX finite element model was initially used to study flexbeam geometries that were 
established by engineering judgment. Because of the indeterminate nature of the problem and 
the millions of possible geometry variations which could be developed, engineering judgment was 
deemed inadequate to quickly define an optimum (lowest stresses) configuration. In addition, 
an optimum configuration for one specific flight condition, such as cruise flight, was significantly 
different from a configuration which would produce minimum stresses in a maneuver condition. To 
address this problem an optimization routine called ADS (Automated Design Synthesis), Ref. (3), 
was added to the HUBFLEX analysis program. The routine permitted holding certain parameters 
of the design fixed or within specified bounds while freeing other variables to achieve a minimum 
stress distribution along the length of the beam. A further enhancement was added by providing 
for the weighted optimization of several flight conditions simultaneously. A block diagram of the 
optimization procedure is provided in Figure (2). 
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Iii practice, the optimization program was used to minimize the maximum stresses (normal and 
shear) caused by maximum maneuver and endurance motions. Design variables were the dimensions 
of the beam element, with upper and lower bounds established for manufacturing feasibility. The 
required damper motion per degree lag motion was selected as a constraint in the optimization 
process. 

The results of this work led to the highly tailored design configurations established for three different 
flexbeams for three different helicopters. These are, ACH (Advanced Composite Hub), advanced 
HARP and most recently, the MDX*flexbeam. The advanced HARP flexbeam was designed in 
less than a week (the old HARP flexbeam took more than 6 months), while it took only half a 
day to design a flexbeam for the MDX main rotor. The reduction of design time from one week 
for advanced HARP to half a day for the MDX flexbeam was purely due to postprocessing of the 
optimization results. The postprocessing of the advanced HARP flexbeam results was carried out 
by hand and relied upon engineering judgment, while the postproccessing of the MDX flexbeam 
optimization results was done by computer and hence human factors and consequently excessive 
time was removed completely' from the optimization process. 

In general, it was found that the endurance level flight condition dictated the geometry of the 
outer portion of the beam where high cycle shear stresses were most critical. The low-cycle high g 
maneuver condition dictated the beam geometry of the inboard end due to high normal stresses. 

Among the three optimized flexbeams, the advanced HARP is the only one which has been fabri- 
cated and successfully tested in fatigue as well as in the Duits Nederlandse Windtunnel (DNW). 


Optimization Results 


The finite element optimization program HUBFLEX/ ADS was used to optimize the advanced 
HARP flexbeam. The composite material selected for the advanced beam was S2-Glass fiber and 
the matrix was Epoxy. The maximum combined normal and/or shear stresses for the one hour 
flight condition as well as endurance were selected as the objective function. To achieve at least 
3% critical damping, a damper motion of 0.1 inch per degree lag was required. Therefore a lower 
constraint was put on the damper motion. In addition to this active constraint, to insure the 
continuity of fibers, a set of constraints was put on the beam cross sectional area variation over 
the length of the beam such that the area was required to either decrease or remain constant, but 
was not allowed to increase from inboard to outboard. In addition to these constraints, lower and 
upper constraints were imposed on the dimensions of the beam element, based on manufacturing 
requirements. 


McDonnell Douglas (X ) 
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The HUBFLEX/ADS computer program was used to optimize the advanced beam based on the 
above optimization statement. A comparison between the initial design and the optimum design for 
the advanced beam is shown in Table(l). As can be seen in this table, it took more than 6 months to 
design the old flexbeam with unacceptable stresses. However with the HUBFLEX/ADS program, 
an acceptable configuration was obtained in a week, with acceptable stresses. The damper motion 
per degree lag, the hinge offsets in flap and lag directions, and the corresponding first frequencies 
for both the old and the advanced beams are acceptable as shown in the table. The alternating 
angles for which the old and the advanced beams were designed are also shown in the table. 

A three-dimensional view of the advanced beam is shown in Figure (3). As seen in the figure, the 
variations of width and thickness are such that the geometric area of the beam section is always 
decreasing or is constant. This particular requirement insures the continuity of the fibers. 


Table(l) 


PARAMETER 

“OLD” HARP 

NEW HARP 

Time to design 
configuration 

6 mo. 

1 wk. 

Normalized 
max combined 
normal stresses 

1.70 

1.04 

Normalized 
max combined 
shear stresses 

0.96 

1.03 

Damper motion 
per deg. lag 

0.18 

0.16 

^ x 100 

5.3 

4.2 

^ x 100 

24.8 

15.6 

ta 

"n 

1.06 

1.04 

u/ < 

n 

0.61 

0.50 

0 

6.8 

6.8 

r 

2.0 

2.0 

6 

11.5 

11.5 


blade attachment 
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HUBFLEX/NASTRAN Comparison 


In the development stage, the HUBFLEX program was validated against 2-dimensional and 3- 
dimensional equivalent beam NASTRAN models and showed excellent correlation. 

For the advanced HARP program, an additional level of verification was obtained by creating 
a 3-dimensional nonlinear anisotropic solid element NASTRAN model of the detailed flexbeam 
structure. The model was constructed using 3048 8-noded (6ctahedrbh)and 8 6-noded (hexahedron) 
solid elements incorporating over 12,000 degrees of freedom, Figure (4). The snubber damper 
centering bearing was rigidly attached to the flexbeam using constraint equations. The chordwise 
and flapwise damper stiffnesses were represented by elastic springs and attached to a rigid element 
pitchcase. The inboard end of the flexbeam was rigidly constrained. The outboard blade attachment 
location was multipoint constrained to allow loads and the pitchcase attachment to be applied at a 
single centerline grid point. 

The detailed nonlinear model was run for flap, lead-lag and feathering one hour flight motions 
using MSC/NASTRAN’s geometric nonlinear solution 64 (whose algorithm is an extension of the 
differential stiffness approach, Ref (4)). The CF loading was applied in the linear elastic step 
and then combined with the corresponding flap, lead-lag and feathering loading in the following 
differential stiffness step. The combined loading was then taken through one additional nonlinear 
interation for an improved solution. Computation time for each run amounted to 8 hours of CPU 
time running on a VAX 11/785. 

The primary loading condition of interest was lead-lag where the HUBFLEX program approximates 
the chordwise stiffness of the split leg/shear web section of the flexbeam by an equivalent ‘I’ beam 
moment of inertia. 



Figure (4) 
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In order to compare with the HUBFLEX output of spanwise displacements, a set of grid points 
on the centerline of the flexbeam was selected to represent the overall deformation. For the com- 
parison of spanwise normal stress distributions, element stress peak values obtained within the 
corresponding cross section location were used. 

Correlation of chordwise loading indicated that the HUBFLEX chordwise stiffness approximation 
produced a stiffer beam resulting in conservative peak stress values that were approximately 10 
percent higher than those of the detailed solid element model. A comparison of the chordwise 
displacement and spanwise stress distribution is provided in Figures (5) and (6). 

For the trivial flapwise and feathering loading cases, correlation was excellent as expected. Peak 
stress values were typically 3 to 5 percent higher for the HUBFLEX model. Figures (7) and (8) 
demonstrate this correlation for the flapwise loading case. 


HUBFLEX vs 3D Solid Nonlinear NASTRAN 


ChordwiM Deflection Comparison 



Flexbeam Station Onchee) 

□ HUBFLEX + Nonfinear NASTRAN 


Figure (5) 

HUBFLEX vs 3D Solid Nonlinear NASTRAN 
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D HUBFLEX Nonfinear NASTRAN 
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Figure (6) 





HUBFLEX vs 3D Solid Nonlinear NASTRAN 


Flapwiee Deflection Comparison 
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Figure (7) 


HUBFLEX vs 3D Solid Nonlinear NASTRAN 


Flapwise Normal Stress Comparison 
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Figure (8) 
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Conclusions 


1. Pitchcase and snubber damper representations are required in the flexbeam model for proper 
sizing resulting from dynamic requirements. 

2. Optimization is necessary for flexbeam design. It reduces the design iteration time and results 
in an improved design. 

3. Inclusion of multiple flight conditions and their corresponding fatigue allowables is necessary 
for the optimization procedure. 

4. A simplified beam model is adequate for flexbeam sizing. The simplified model gives excel- 
lent accuracy (compared to a detailed 3D nonlinear NASTRAN analysis) with a significant 
reduction in computational time. 

5. HUBFLEX’s rapid computation capability enables design parameters to be easily modified 
and implemented. 
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